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ABSTRACT 

Using data from the Sloan Digital Sky Survey (SDSS; data release 7), we have conducted a search for 
local analogs to the extremely compact, massive, quiescent galaxies that have been identified at z > 2. 
We show that incompleteness is a concern for such compact galaxies, particularly for low redshifts 
(z < 0.05), as a result of the SDSS spectroscopic target selection algorithm. We have identified 63 
M* > lO 16 ' 7 M Q (w 5 x 10 10 M ) red sequence galaxies at 0.066 < z spC c < 0.12 which are smaller than 
the median size-mass relation by a factor of 2 or more. Consistent with expectations from the virial 
theorem, the median offset from the mass-velocity dispersion relation for these galaxies is 0.12 dex. 
We do not, however, find any galaxies with sizes and masses comparable to those observed at z ~ 2.3, 
implying a decrease in the comoving number density of these galaxies (at fixed size and mass) by a 
factor of > 5000. This result cannot be explained by incompleteness: in the 0.066 < z < 0.12 interval, 
we estimate that the SDSS spectroscopic sample should typically be > 75% complete for galaxies with 
the sizes and masses seen at high redshift, although for the very smallest galaxies it may be as low 
as ~ 20%. In order to confirm that the absence of such compact massive galaxies in SDSS is not 
produced by spectroscopic selection effects, wc have also looked for such galaxies in the basic SDSS 
photometric catalog, using photometric redshifts. While we do find signs of a bias against massive, 
compact galaxies, this analysis suggests that the SDSS spectroscopic sample is missing at most a few 
objects in the regime we consider. Accepting the high redshift results, it is clear that massive galaxies 
must undergo significant structural evolution over z < 2 in order to match the population seen in the 
local universe. Our results suggest that a highly stochastic mechanism like major mergers cannot be 
the primary driver of this strong size evolution. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: fundamental parameters 



1. INTRODUCTION 

In the simplest possible terms, the naive expectation 
from hierarchical structure formation scenarios is that 
the most massive galaxies form last. This is in contrast 
to the observation that the bulk of cosmic star formation 
occurs in galaxies wit h progressively lower s tellar masses 
at later times (e.g. Uuneau et all 120051 : IZheng et al.1 
l2007t iDamen et al.ll2008[) : the so-called downsizing of 
galaxy growth. These observations have been accommo- 
dated within the ACDM framework with the introduc- 
tion of a quenching mechanism (e.g. iMenci et alj [20051 : 
ICroton et all l2006UCattaneo et all T2008f) . which oper- 
ates to shut down star formation in the most massive 
galaxies; this mechanism is also required to correctly 
pre dict the absolute and re l ative numbers of red galax- 
ies (iDekel fc Birnboimll2006l:lBerl et aIll2007l : lFaber et all 
l2007h . With this inclusion, models thus predict that a 
significant fraction of the local massive galaxy popula- 
tion should have finished their star formation relatively 
early in the history of the universe, with later mergers 
working to build up the most massive galaxies. 

There is thus a crucial distinction to be made between a 
galaxy's mean stellar age, and the time si nce that galaxy 
has a ssumed its present form (see, e.g., IDe Lucia et al.l 
2006) : the most massive galaxies are expected to be both 
the oldest and the youngest galaxies. They are the oldest 
in the sense that their progenitors are expected to form 



first in the highest cosmic overdensities — however, these 
stars are only assembled into their z = configuration 
relatively recently. 

This leaves (at least) two open questions relating to 
the quenching of star formation and the formation and 
evolution of massive galaxies: 1.) When does star forma- 
tion stop in massive galaxies, and 2.) What happens to 
galaxies after they have stopped forming stars? 

In connection with the first of these questions, deep 
spectroscopic surveys have identified massive galaxies 
with little or no ongoing star formation at 1 < g 5, 2 
(e.g. ICimattil 120041: (Glazebrook et all 12004 IMcCarthvl 
2004a; Daddi et al.ll2004l ). At th e same time, colo r se- 
lection techniques like the ERO (IMcCarthvl I2004E and 
references therein) , DRG ([Franx et alj |2003T ). or BzK 
(|Daddi et alj 12005) criteria have been used to identify 
massive, passive galaxies at high redshifts. While these 
techniques are deliberately biased towards certain kinds 
of galaxies and certain redshift intervals, advances in 
techniques for photometric redshift estimation and stel- 
lar population modeling have allowed the selection of 
mass-limited samples, and so the construction of rep- 
resentative sam ples of the high redshift massive galaxy 
population (e.g. Ivan Dokkum et aTll2006h . 

By obtaining very deep rest-frame optical spectra of 
a photometri c-redshift se l ected s ample of massive galax- 
ies at z > 2, iKriek et al.l (|2008af ) made a significant ad- 
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vance on previous spectroscopic and photometric stud- 
ies. Of the 3 6 z S pec > 2, M* > 10 11 M galaxies in the 
iKriek et al] (|2008al) sample, 16 were shown unambigu- 
ously to have evolved stellar populations and little or no 
ongoing star formation. These galaxies also seem to form 
a red sequen ce in (B — V) color although at low signif- 
icance f3.3cr; iRriek et al.1 [20 08b). In other words, these 
massive galaxies appear both to have assembled stellar 
populations similar to galaxies of comparable mass in 
the local universe, and to have had their star formation 
effectively quenched. 

Using Keck laser guide-star assisted adaptive optics 
and Hubble Space Telescope imaging, van Dokkum et 
al. (2008, hereafter vD08) measured sizes for 9 of the 16 
strongly quenched galaxies from the IKriek et al.l (|2008af) 
sample. They found (rest-frame optical) effective radii 
in the range 0.5 — 2.4 kpc; that is, smaller than typi- 
cal galaxies of the same mass in the local universe by 
factors of 3 — 10. These galaxies have stellar mass den- 
sities, measured within the central 1 kpc, which are 2 — 
3 tim es higher than typica l l ocal galaxies of the same 
mass (jBezanson et al.ll2009H . ICimatti etldl ()2008[ ) and 
Damjanov et al. (2009, hereafter D09) have found sim- 
ilarly compact sizes for massive galaxy samples drawn 
from 1 < z < 2 spectro s copic surveys. Further, 
Ivan Dokkum. Kriek fe Franx! (|2009h have recently mea- 
sured a velocity dispersion of SlO^gg 5 km/s for one of the 
galaxies in the vD08 sample, based on a 29 hr NIR spec- 
trum; this extremely high value is consistent with the 
galaxy's measured mass and size. (See also Cappellari et 
al. 2009.) 

By providing rest-frame optical size measurements 
for a representative, mass-limited sample of galaxies 
spectroscopically-confirmed to have little or no ongoing 
star formation and z > 2, these re s ults confirm and con- 
solida te the work of [Daddi et all (12005D. ITruiillo et all 
(I2006h. iTruiillo et al.l (|2007f l. IZirm et alJ (|2007t l. and 
iToft et all (l2007f). as well a s 1 < z < 2 result s from , 
e.g.. lLonghetti et all (120071) and ISaracco et all (I2009D, 
and z < 1 results from Ivan der Wei et al.l ( 2008f) . (See 



also Buitrago et al., 2008.) 

The significance of these results is that while the mas- 
sive and largely quiescent galaxies at z > 2 have stellar 
populations that are consistent with their being more or 
less 'fully formed' early type galaxies, they must each 
undergo significant structural evolution in order to de- 
velop into galaxies like the ones seen in the local uni- 
verse. Taken together, these results thus paint a con- 
sistent picture of strong size evolution among massive, 
early type and/or red sequence galaxies 1 — both as a pop- 
ulation and individually — ev en after their star formation 
has been quenched (see also lFranx et al.l r2008f ). What- 
ever the mechanism for this growth in size (see, e.g. , 
Fan et al.l 120081: iHopkins et all l2009t iNaab et alJ 120091: 
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galaxies is not monolithic. 

The aim of this paper is to test the proposition that 
there are no galaxies in the local universe with sizes and 

1 There is considerable, but not total, overlap between color- 
selected samples of red sequence galaxies, and morphology-selected 
samples of early type galaxies. While it is common to use these 
terms as if they were more or less interchangeable, it should be 
remembered that they are not. 



masses comparable to those found at z > 1.5 — this is the 
crux of the argument against the monolithic formation 
of massive galaxies. This work is based on the latest 
data products from the Sloan Digital Sky Survey (SDSS; 
lYork et alJfeOOOt IStrauss et all 12002). In particular, we 
will focus on the possibility that such galaxies have been 
overlooked in SDSS due to selection effects associated 
with the construction of the spectroscopic target sample. 

The structure of this paper is as follows: We describe 
the basic SDSS data that we have used in SJ21 In S3 
we define our sample of compact galaxy candidates, and 
present several checks to confirm that these galaxies are 
indeed unusually small for their stellar masses. Then, 
in 21 we consider the importance of the SDSS spectro- 
scopic selection for massive, compact galaxies. In this 
Section, we also compare ou r z ~ 0.1 compact galaxy 
candidates with the vD08 and lD09l samples. In Appendix 
IA1 we provide a complementary analysis in order to con- 
firm our conclusion that the apparent differences between 
the high- and low-redshift samples cannot be explained 
by selection effects, including an estimate for the num- 
ber of compact galaxies that may be missing from the 
SDSS spectroscopic sample. Finall y, in j}5l we compare 
our results to a similar studi es by iTruiiilo "eTall (|2009h 
and IValentinuzzi et aT] (|2009t ) , and briefly examine the 
properties of our compact galaxies' stellar populations in 
comparison to the general z ~ 0.1 red sequence galaxy 
population. A summary of our main results is given in 
S|6] Throughout this work, we assume the concordance 
cosmology; viz.: fl m — 0.3, fl\ = 0.7, and H n = 70 
km/s/Mpc. 

2. BASIC DATA AND ANALYSIS 

The present work i s based on Data Release 7 (DR7; 
Abazajian et~aT][2009] ) of the SPSS, accesse d via the Cat- 
alog Archive Server 2 (CAS; iThakar et aTll2008h . In this 
section, we describe the basic SDSS data that we have 
used, and our analysis of it. We will search for compact 
galaxy candidates in the SDSS spectroscopic catalog; to 
this end, we will only consider sciencePrimary objects 
(a flag indicating a 'science-grade spectrum, and weeding 
out multiple observations of individual objects) with ei- 
ther a 'star' or 'galaxy' photometric type (ie., a genuine 
astronomical source). The details of the SDSS spectro - 
scopic sample selection are given in IStrauss et al.1 (|2002h : 
we will summarize the most relevant aspects of this pro- 
cess in i)4.1l 

2.1. The Basic SDSS Catalog 

For the basic SDSS catalog, there are two different 
methods for performing photometry. The first, the 'Pet- 
rosian' magnitude, is derived from the observed, az- 
imuthally averaged (ID) light profile. The Petrosian 
radius is defined as the point where the mean surface 
brightness in an annulus drops to a set fraction (viz. 0.2) 
of the mean surface brightness within a circular aper- 
ture of the same radius. Within SDSS, the Petrosian 
aperture is defined to be twice the Petrosian radius; this 
aperture will contain 99 % of the total light for a well re- 
solved exponential disk, but may miss as much as 18 % of 
the light for a de Vaucoul eurs i? 1 / 4 profile (|Strauss et al.1 
l200l iBlanton et al.ll200l . 
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The second photometric measure is derived from fits to 
the observed (2D) distribution of light in each band, us- 
ing a sector-fitting technique, in which concentric annuli 
are di vided into 12 3 ° sect ors, as described in Appendix 
A.l of iStrauss et al.l (j2002j ). These fits are done assum- 
ing either an exponential or a de Vaucouleurs profile, 
convolved with a fit to the appropriate PSF. For each 
profile, the structural parameters (viz. axis ratio, posi- 
tion angle, and scalelength) are determined from the r 
band image. The more likely (in a \ 2 sense) of the two 
profile fits is used to define 'model' magnitudes for each 
galaxy. For the ugiz bands, these parameters are then 
held fixed, and only overall normalization (ie. total flux) 
is fit for. 

The basic catalog also provides two different measures 
of size, associated with these two magnitude measure- 
ments. The Petrosian half-light radius, R$q, is defined as 
the radius enclosing half the 'total' Petrosian flux. The 
catalog also contains best fit structural parameters, in- 
cluding the effective radius, from a separate set of fits to 
each band independently, again for both an exponential 
and a de Vaucouleurs profile. Note that whereas the Pet- 
rosian magnitude and size are derived from the observed, 
PSF-convolved radial profile, the model values provide a 
PSF-corrected measure of the intrinsic size. 

We use model magnitudes to construct ugriz SEDs 
for each object, since these measurements are seeing- 
corrected. From DR7, the basic SDSS photometric cali- 
bration has been refined so that the photometry is given 
in the AB magnit ude system without the ne ed for any 
further corrections (iPadmanabhan et al. 2008) . For mea- 
suring sizes, we will rely on the best-fit model effective 
radius, R c , as determined from the z band. We also 
adopt a minimum measured size of 0'.'75, corresponding 
to half the median PSF FWHM for the SDSS imaging; 
we will plot all galaxies with observed sizes smaller than 
0'.'75 as upper limits. (None of our conclusions depend 
on the choice of this limit, which ultimately affects only 
5 of our lower-mass compact galaxy candidates.) 

2.2. Derived Quantities 

We have derived rest-frame photometry for each ob- 
ject, based on its observed ug riz SEP and redsh ift, us- 
ing the IDL utility InterRest (jTavlor et al.ll2009f ), with 
a redshift grid of Az = 0.001. In order to minimize 
the k-corrections and their associated errors, we deter- 
mine rest-frame photometry through the ugriz filters 
rcdshifted to z — 0.1, which we denote with a super- 
script 0.1. We estimate that the systematic uncertainties 
are at the level of < 0.02 mag. The agreement between 
our interpolated rest-frame photomet ry and that derived 
using the SDSS kcorrect algorithm (jBlanton fc Roweisl 
2007) is very good: our derived (u — r) colors are ~ 0.02 
mag bluer for blue galaxies, and ~ 0.03 mag redder for 
red galaxies. 

We make use of stellar mass estimates provided by the 
MPIA Garching group. 3 JB has fit the ugriz photom- 
etry of all galaxies using the synthetic s tellar popula- 
tio n library described by IGallazzi et alj (|2005l ). based 
on iBruzual fc CharloU (2003) models and assuming a 
Chabrier ( 2003) stellar init ial mass function (IMF ) in the 
range 0.1—100 M Q . The IGallazzi et al.l (|2005ft library 

3 Available via http://www.mpa-garching.mpg.de/SDSS/DR7 . 



contains a large number of Monte Carlo realizations of 
star formation histories, parameterized by a formation 
time (1.5 < t{ olm /[Gyr] < 13.5), an exponential decay 
rate (0 < 7/[Gyr~ 1 ] < 1), and including a number of 
random star formation bursts (randomly distributed be- 
tween tf orm and 0, normalized such that 10 % of galaxies 
experience a burst in the last 2 Gyr). In the fitting, the 
photometry has been corrected for emission lines under 
the assumption that the global emission line contribution 
is the same as in the spectroscopic fiber aperture. 

The agree ment between thes e SEP- fi t mass estimates 
and those of iKauffmann et all (|2003a[ L which were de- 
rived from spectral line indices, are excellent: the median 
offset is -0.01 dex, with a scatter on the order of 0.1 dex. 
For the highest masses, however, the SEP-fit results are 
slightly less robust: for M* > 10 11 M©, the median for- 
mal error is < 0-10 dex, compared to < 0.06 dex for the 
IKauffmann et all (|2003af ) estimates. 

In the upper panel of Figure[TJ we show the stellar mass 
to light ratios, M*/Li, for 0.066 < z < 0.12 galaxies as 
a function of their 01 (g — i) color; here again, Li should 
be understood as referring to the i-band filter redshifted 
to z — 0.1, or 01 i. Notice that, at least for these mass 
estimates, M*/L is a very tight function of color. In the 
main panel of this Figure, the red line shows the median 
M*/Li in narrow color bins. Making a simple linear fit 
to these points, we find: 

]og(M*/Li) = -0.82 + 0.83 x °- 1 (g - i) , (1) 

where both M* and Li are in solar units. (The absolute 
magnitude of the sun in the i band is 4.58.) This re- 
lation is shown in Figure [1] as the solid blue line. We 
present this relation as an alternative to the popular 
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Fig. 1.— The mass-to-light ratios, M,/Li, of 0.066 < z < 0.12 
galaxies as a function of their (g — i) color. (Here, M t /L{ is under- 
stood to relate to the i band redshifted to z = 0.1.) The greyscale 
shows the (linear) data density in cells, where the data density is 
high. In the main panel, the red line shows the median M*/Li 
in narrow color bins; the blue line is a linear fit to these points. 
In the lower panel, we have simply subtracted away the median 
relation; in this panel, the error bars show the 16/84 percentiles in 
color bins. The simple linear relation shown provides an acceptable 
description of the observed relation, to within 0.02 — 0.04 dex; the 
global RMS offset from this relation is 0.032 dex. In order to avoid 
selecting 'catastrophic failures' in terms of stellar mass estimates, 
we will consider only those galaxies that fall within 0.25 dex of the 
median M./L;- - 1 (g - i) relation, and with 0.4 < 0,1 (g- i) < 1.8, 
as shown by the box in the lower panel. 
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Fig. 2. — The size-mass relation for massive, red sequence galaxies showing the importance of the SDSS spectroscopic selection criteria. 
The points show SDSS galaxies (0.066 < z < 0.12) selected to have 01 (u — r) > 2.5. The yellow points show the median size in narrow 
bins of stellar mass; t he error bars show the 16/84 percentiles in each bin. A fit to this median size— mass relation for red sequence galaxies 
is consistent with the Shcn ct al. (2003) relation for early type galaxies (dashed line), albeit offset by 0.05 dex. Individual galaxies that 
we have visually inspected (M* > 10 10 ' 7 Mq; Alogi? e < —0.3 dex) are marked with large symbols. Galaxies with M/Ls that differ 
significantly from the main color— M/L relation shown in Figure [T] are marked with small blue crosses. Galaxies with obvious problems 
in their photometry (especially those affected by the presence of a bright nearby star or blended with other galaxies) are marked with a 
small red cross; those that look okay are plotted as circles. Galaxies with observed sizes smaller than 0'.'75 are plotted as upper limits, 
assuming a size of 0"75. The different lines show how the principal selection limits for spectroscopic followup translate onto the (M,, R e ) 
plane for z = 0.12, 0.10, 0.066, 0.050, and 0.35 (top to bottom): the diagonal, long-dashed lines show the star/galaxy discrimi nator ; the 
short-dashed boxes show the 'saturation' selection limit, and the diagonal dotted lines show the' cross-talk' selection limit (see i!4. II for a 
detailed discussion). Galaxies lying below these lines will not be targeted spectroscopically. 



iBell fc de Jonel (|200lf) or lBell et all ((20031 ) relations. 

In the lower panel of Figure [TJ we show the disper- 
sion around the median relation; in this Figure, the error 
bars show the 16/84 percentiles in bins of color. Over- 
all, the dispersion around this relation is just 0.032 dex. 
Note that while the simple linear relation given above 
provides an acceptable description of the 'true' relation, 
systematic offsets exist at the 0.02 — 0.04 dex level. The 
global mean and random offset from this linear relation 
are 0.002 dex and 0.040 dex, respectively. 

In both panels, the small grey pluses show points that 
fall outside the plotted range. Notice that there are 
a small fraction of galaxies that fall well off the main 
M*/L-color relation, some by an order of magnitude or 
more. These galaxies also lie significantly off the main 



stellar mass-dynamical mass relation and are very likely 
to represent catastrophic failures of the stellar mass SED- 
fitting algorithm. This presents a problem when it comes 
to looking for outliers in the mass-size plot: selecting 
the most extreme objects may well include those ob- 
jects with the largest errors. For this reason, we will 
restrict our attention to those objects that fall within 
0.25 dex (« 7.8a) of the main M*/L-color relation, and 
with 0.2 < 0A (g — i) < 1.8, as shown by the box in 
the lower panel of Figure [TJ This selection excludes just 
under 600 of the 223292 galaxies shown in Figure [T] 

3. SEARCHING FOR MASSIVE, COMPACT, 
EARLY-TYPE GALAXIES IN THE LOCAL 
UNIVERSE 
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Fig. 3. — Illustrative examples of the galaxies we consider. Clockwise from the top-right, we show a 'normal', massive early type galaxy 
that lies very close to the median size-mass and velocity dispersion-mass relations, two compact galaxy candidates where visual inspection 
suggested problematic size measurements, and three of our compact galaxy candidates. Each of the three compact galaxy ca ndid ates shown 
in this Figure have observed velocity dispersions that are approximately consistent with their small measured sizes (see ^ 3- 3 D . For each 
galaxy, wc show the SDSS SkyServer thumbnail image used for visual inspection, as well as the galaxies' observed spectra; the boxes show 
the SEDs from the photometry, scaled to match the spectroscopic flux in the r band. 



3.1. Identifying Massive, Compact Galaxy Candidates 

Figure [2] shows the size-mass plot for a sample of mas- 
sive, red-sequence galaxies drawn from the SDSS DR7 
spectroscopic sample. These galaxies have been selected 
to have 0A (u - r) > 2.5 and 0.066 < z < 0.12. These 
rcdshift limits have been chosen to minimize the impor- 
tance of selection effects and measurement biases, which 
we will discuss in ij4.ll For now, we note that, mapping 
the r < 17.77 spectroscopic limit onto M*(z), we should 
be highly complete (volume limited) for M* > 10 10 7 
Mq and z < 0.12. As a very simple check on this, we 
note that for this sample, the median redshift in narrow 
bins of stellar mass is within the range z — 0.098 — 0.102 
for all M* > 10 10 7 M Q ; the volumetric center of the 
0.066 < z < 0.12 bin is z = 0.10. 

The yellow points in this Figure show the median size 
in narrow bins of stellar mass; the error bars show the 
14/86 percentiles. For comparison, the long-dashed line 
show s the local s i ze-ma ss relation for early-type galaxies 
from IShen et al.l ((2003), corrected for differences in as- 
sumed IMF and cosmolo gy. Contrary to the findings of 
IValentinuzzi et ail (|2009t ). a simple fit to the size-mass 
relation for red sequence galaxies ( 01 (it — r) > 2 .5) sh own 
in Figure [5] is consistent with the IShen et ail (|2003| ) re- 
lation for early type (n > 2.5) galaxies, albeit offset in 
size by 0.05 dex or, equivalently, by —0.09 dex in mass. 
At fixed mass, the mode of the distribution is similarly 
offset (see Figure [7]); this does not appear to be due to 
large numbers of late type galaxies in the sample. 



We next select and study very compact galaxies from 
within the red sequence sample shown in Figure O At 
first glance, it appears that there may be a few galaxies 
that lie well below the main size-mass relation. How- 
ever, it must be remembered that by selecting the most 
extreme outliers, we will also be selecting those objects 
with most egregious measurement errors. 

For this reason, we have individually visually inspected 
all M» > 10 10 7 M Q galaxies with inferre d sizes that 
are les s than half the size predicted from the IShen et all 
(2003) relation; it. AR C < -0.3 dex. For sizes smaller 
tha n the median relat ion, the distribution of sizes around 
the IShen et all (|2003t ) relation is very well described by 
a Gaussian with a = 0.11 dex; this AR C cut thus corre- 
sponds to selecting those galaxies whose sizes are smaller 
than the mean size (at fixed mass) at the > 2.7cr level. 
(Adopting our own fit to the size-mass relation, this se- 
lection translates to AR C < —0.35 dex; our results are 
otherwise unchanged.) 

We have inspected 280 such objects, and discarded 
those where there are obvious reasons to distrust the size 
measurements. The most common reasons for discarding 
galaxies were confusion with other galaxies (99 galaxies, 
including 19 good merger candidates, and two possible 
lenses), or with the extended halos, diffraction spikes, 
and/or reflections of bright stars (62 galaxies). A fur- 
ther 19 galaxies were clearly disk-like, 5 showed marked 
asymmetries, and 1 had a very strong AGN spectrum; 
these candidates were also discarded. We discarded a 
further 3 objects with bad or missing data. 
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TABLE 1 

Properties of our Compact Galaxy Candidates 



RA 


dec 


z 


(u — r) t, s 


O-^u-r) 


logM, 


650 z 


71 


©50 z 


Re 


cr 


T 


Z 


(1)« 


(2) a 


(3) a 


(4)o 


(5) 


(6) 6 


(7)° 


(8) c 




(10) 


(ll) 6 


(12) d 


(13) d 


190.16652 


13.81563 


0.0855 


2.577 


2.642 


10.700 


0.912 


4.31 


0.917 


1.462 


160 






127.02722 


55.37988 


0.0665 


2.828 


2.999 


10.701 


1.152 


3.26 


1.211 


1.468 


191 






225.31708 


30.58266 


0.0980 


2.858 


2.877 


10.705 


0.808 


3.87 


0.773 


1.464 


195 






227.08531 


7.25325 


0.0764 


2.963 


3.113 


10.709 


0.929 


5.13 


1.188 


1.345 


199 






215.41043 


40.03233 


0.0998 


2.803 


2.813 


10.709 


0.795 


4.11 


0.670 


1.464 


176 


9.775 


0.035 


222.12988 


26.48791 


0.1058 


2.549 


2.540 


10.712 


0.750 


4.65 


0.722 


1.453 


155 






118.81702 


33.22864 


0.0980 


2.680 


2.697 


10.713 


0.803 


2.67 


0.739 


1.454 


154 


-99 


-99 


143.05707 


11.70454 


0.0811 


2.690 


2.776 


10.726 


0.750 


5.52 


0.830 


1.146 


166 


9.255 


0.132 


204.66577 


59.81854 


0.0704 


2.857 


3.012 


10.731 


0.943 


3.92 


0.974 


1.267 


235 


9.845 


0.229 


230.28553 


24.21978 


0.0809 


2.922 


3.028 


10.733 


0.952 


5.90 


1.298 


1.453 


153 







Note. — We give only the first ten candidates here; the properies of the full sample of 63 galaxies is available as an electronic table via 
|http://www.strw.leidenuniv.nl/~ent/| Col.s (1) and (2): postition in decimal degrees (J2000); Col. (3) spectroscopic redshift; Col.s (4) 
and (5): observed and rest-trame colors; Col. (6) stellar mass in units of solar masses; Col. (7) apparent De Vaucouleurs effective radius in 
arcsec; Col. (8) and (9) Sersic index and apparent Sersic effective radius in arcsec; Col. (10) physical De Vaucouleurs effective radius in kpc; 
Col. (11) velocity dispersion in km/s; Col. (12) luminosity-weighted age in Gyr; Col. (13) mean metallicity. In Col.s (12) and (13), data is 
only given for those objects that appear in DR4; values of -99 indicate unsuccessful fits to the spectra. SOURCES. — a : the default SDSS 
llYork et al.ll2000l ; [StFauss et al.ll2002h catalog for DR7 HAbazaiian et al.|[2009l) . accessed via CAS l IThakar et al.|[2008T) ; b : the MPA-JHU 
catalog for DR7 (accessible via http: //www.mpa-garching .mpg.de/SDSS/DR7/); c : the NYU VAGC l IBlanton et al.1120051) for DR7; d : the 
stellar age and metallicity catalog ot Gallazzi et al. (2005) tor JJR4. 



In Figure 02 we show several illustrative examples of 
the galaxies we are considering. On the right-hand side 
of this Figure, we show a 'normal' early type galaxy, with 
M» ~ 10 11 Mq, which falls very close to the IShen et al.l 
(2003) relation. Below this, we show two of the compact 
galaxy candidates that we have rejected on the basis of 
visual inspection. On the left-hand side of this Figure, we 
show three of the compact galaxy candidates of different 
stellar masses that we have retained after visual inspec- 
tion. For each galaxy, we show the thumbnail image from 
the SDSS SkyServer 4 , used for visual inspection. We also 
show each galaxy's observed spectrum and photometry; 
here, we have scaled the photometry to match to the 
integrated r-band flux from the observed spectrum. 

In addition to these galaxies with suspect size mea- 
surements, we have excluded a further 27 galaxies whose 
SED-fit M*/Ls are offset from the main color-M*/_L re- 
lation shown in Figure [T] by more than 0.25 dex. If we 
use Equation [1] to derive new stellar mass estimates for 
these galaxies, all of these galaxies move back into the 
main cloud in both Figure[2]and a stellar mass-dynamical 
mass plot, with mean/median offsets of < 0.02 dex in 
both cases. 

The 190 galaxies discarded on the basis of inspection 
are shown in Figure [2] as small red crosses; the small blue 
crosses show the 27 galaxies with discrepant M*/Ls. As 
a function of Ai? e , the fraction of inspected sources that 
have been discarded goes fairly smoothly from 60 % for 

AR C 0.3 dex to ~ 100 % for AR C < -0.5 dex. The 

discarded fraction has a similar dependence on mass: it 
is - 60 % for M* ~ 10 10 7 M Q , rising to ~ 85 % for 
M* ~ 10 11 Mq, and 100 % for M* > 10 114 M . 

This leaves us with a sample of 63 massive, compact, 
early- type and red sequence galaxy candidates; these are 

4 Also accessible via CAS at |http: / / c as.sdss.org 



are marked in Figure[2]with heavy black circles. Of those 
galaxies that we have retained, 8 % (5/63) have observed 
sizes smaller than 0'.'75; all of these have < 10 11 M . 
We have provided the properties of our compact galaxy 
candidates in Table [TJ 

3.2. Are the Size Measurements Wrong? 

We have performed a number of checks to validate the 
small measured sizes of our compact galaxy candidates. 
The compact galaxy candidates do not have significantly 
larger size measurement errors in comparison to the full 
sample shown in Figure [2] For both the r- and z-bands, 
our candidates are not anomalous in a plot of Petrosian 
half-light radius versus model effective radius, nor are 
they anomalous in a plot of r-band size versus z-band 
size. For all but two of the candidates, the Petrosian 
and model magnitudes agree to within 0.15 mag. The 
mean offset between model and Petrosian magnitudes is 
-0.06 mag for the compact galaxies, compared to -0.08 
mag for the full sample shown in Figure [2] That is, the 
compact candidates appear to be well described by the 
de Vaucouleurs model fits. 

For the New York Un iversity (NYU ) Valu e Added 
Galaxy Catalog (VAGC), iBlanton et all (gll) have fit 
the radially-averaged light profiles of each object, fitting 
for the Sersic index as a free parameter over the range 
< n < 6. In order to explore further the issue of the 
quality of the de Vaucouleurs profile fits, we have gone 
to the NYU VAGC for DR7, and looked up the Sersic fit 
results for each of our candidates. 

In Figure [H we show the distribution of Sersic param- 
eters for our candidates, as well as a comparison between 
the Sersic and de Vaucouleurs sizes. First, we note that 
nearly all (59/63) of our compact galaxy candidates have 
n > 3; these are not late- type (exponential) galaxies. It 
is therefore unsurprising — but still reassuring — that the 
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Fig. 4. — Comparison between effective radii derived assuming 
a de Vaucouleurs (n = 4) profile and assuming a Sersic profile 
(0 < re < 6). Whereas the basic SDSS catalog uses a sector- fitting 
technique to fit either an exponentia l (re = 1) or a de Vau couleurs 
(re = 4) profile, for the NYU VAGC, [Blanton et al.l (120051') have fit 
the radial profiles of each object assuming a general Sersic model 
(0 < re < 6). This Figure shows the ratio of these two sizes for our 
compact galaxy candidates, based on the z-band data, as a function 
of Sersic index re. Almost all candidates have re > 3 — these galaxies 
are not obviously exponentials. However, ap proximately 25 % have 
re = 5.9; the maximum value allowed in the Blanton et al. (2005) 
fits. For these galaxies, the median ratio between the two size 
measurements is 0.88, with an RMS scatter of 0.1 dcx. In the 
main panel, we show a least-squares fit to the data, the dispersion 
around this relation is < 0.1 dex. 



two size measures agree quite well: for the median galaxy 
among our candidates, the de Vaucouleurs size is ~ 10% 
smaller than the Sersic size; the RMS dispersion is 0.10 
dex. For comparison, the median quoted error for the de 
Vaucouleurs size measurements is 4.6%. 

Notice that about a quarter (17/63) of our candidates 
have n = 5.9 in the NYU VAGC; this is the maximum 
value allowed in the fits. These galaxies are consider- 
ably more centrally-concentrated than the canonical de 
Vaucouleurs i? 1 / 4 -law profile. However, the trend with 
increasing Sersic index is for the de Vaucouleurs size, 
RueVi to be systematically lower than the Sersic size, 
i?Scr: making a least-squares fit to the data shown in 
Figure H we find log -R De v/#Ser = -0.02 - 0.05 (n - 4). 
If these galaxies do have n > 6, then we may well be 
u nderestima t ing th eir sizes by > 25 %. 

IGuo et al . (2009) have recently demonstrated that as 
a result of biases in the way the background sky level 
is estimated for the Sersic fitting, the NYU- VAGC sizes 
are systematically underestimated at the > 15 % level 
for n > 5. This problem is progressively worse for large 
sizes (9 C > 1") and bright magnitudes (r < 16); for 
our compact galaxy candidates, the effect is likely to be 
at the ~ 20 % level. But note this implies that the 
difference between the de Vaucoleurs and Sersic sizes is 
even greater than Figure [4] might suggest: the sizes of 
the n > 5 compact galaxies may be underestimated by 
as much as > 30 %. 

As a final check, therefore, we have also re-derived 



Sersic effective r adii for our comp act galaxy candidates 
using GALFIT (jPeng et al.l 120021 ) and done a similar 
comparison as for the NYU VAGC sizes. The agree- 
ment between the GALFIT and VAGC Sersic indices is 
quite good, with an rms difference in n of 1.1. Again 
the vast majority of objects have n > 3. There are 19 
objects that are assigned the maximum allowed value of 
n = 8, but only 9 of these have n — 5.9 in the VAGC. 
Making a similar fit to the difference between the default 
De Vaucouleurs and the GALFIT Sersic effective radii, 
we find logiJDev/fisor = 0.08 - 0.08 (n - 4). As before, 
we may be underestimating the sizes of high n galaxies 
by 10 — 35 %, although this comparison suggests that we 
may also be overestimating the sizes of the few candi- 
dates with n < 4. The median galaxy has a GALFIT 
Sersic effective radius 15 % smaller than the default De 
Vaucouleurs value. Lastly, we note that there is a definite 
mass-dependence to the agreement between the GALFIT 
Sersic and default De Vaucouleurs effective radii, such 
that all but one of the galaxies for which the sizes agree 
to within 20 % have M* > 10 11 M . 

To summarize the results of this section, then: com- 
parison to ID and 2D Sersic fits does not suggest that the 
default De Vaucouleurs effective radii from the SDSS cat- 
alog are catastrophically wrong for our compact galaxy 
candidates; if anything, these comparisons suggest that 
we may in fact be underestimating the sizes of these 
galaxies by 10 — 30 %. 

3.3. A Consistency Check Based on Velocity 
Dispersions 

Assuming that elliptical galaxies are structurally self- 
similar, the virial theorem implies that M* cx R e a 2 . At 
fixed mass, galaxies with small sizes should therefore 
have higher velocity dispersions, with Aa cx (Ai? c ) -1 / 2 . 

In order to verify that the observed velocity disper- 
sions of our compact galaxy candidates are consistent 
with their being genuinely small, in the lower panel of 
Figure [5] we plot the offset from the local size-mass rela- 
tion for early type galaxies, Alogi? G , against the offset 
from the M*-cr relation, Alogcr; the M*-cr relation it- 
self is shown in the upper panel of the Figur e. For the 
lower panel of this plot, we have shifted the IShen et al.l 
(2003) relation upwards in size by 0.05 dex to be con- 
sistent with the present data set; our conclusions do not 
depend on this decision. The greyscale and points show 
those 0.066 < z < 0.12 galaxies with °- 1 (u — r) > 2.5 and 
ill* > 10 10 ' 7 M Q ; the red circles indicate those galaxies 
that we have identified as compact. 

For our compact galaxy candidates, the median offset 
from the size-mass relation is Alogi? = —0.38 dex. We 
would therefore expect a median offset from the M*-cr 
relation relation of Alogcr = —0.5 x —0.38 = 0.19 dex. 
The median value for Alogcr is 0.12 dex — roughly 85 % 
of the expected value, and ~ 1.5 times greater than the 
intrinsic scatter in the relation. Overall, these results 
are fairly consistent, although they do indicate that the 
sizes may be underestimated and/or the masses may be 
overestimated at the level of 10-20 %. We note that 
the difference between the default SDSS and the NYU 
VAGC size measurements can account for at least half of 
this effect (see E f3~2]) . 

There is one of our compact galaxy candidates how- 
ever, whose velocity dispersion is clearly inconsistent 
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Fig. 5. — Using observed velocity dispersions to validate the 
measured sizes of our compact galaxy candidates. Upper panel: 
the mass-velocity dispersion relation for red sequence galaxies at 
0.066 < z < 0.12. The points and greyscale show the SDSS data. 
The yellow plusses show the median velocity dispersion in narrow 
bins of stellar mass; the solid yellow line shows a simple fit to 
these points for 10.7 < logM* < 11.5. The red circles highlight 
our compact galaxy candidates. Lower panel: the offset from the 
M*— R e relation, plotted against the offset from the M»-cr relation 
for M* > 10 10 - 7 Mq galaxies with - r) > 2.5. If the offsets 
from these two relations is a function of galaxy size, then we ex- 
pect A logo- = —0.5 X A log R e (long dashed line). Our compact 
galaxy candidates are shown as the red circles. In general, the 
observed velocity dispersions support the idea that our compact 
galaxy candidates are indeed relatively small; this is particularly 
true for those with a > 200 km/s. There is one clear exception, 
marked with a cross in both panels; this galaxy is also the most 
extreme outlier in Figure [4] 



with its being massive and compact, which we have 
marked in Figure [5] with a cross; indeed, it has the low- 
est observed velocity dispersions of all of our compact 
galaxy candidates. This galaxy is also the biggest outlier 
in Figure |U We will discuss this object in more detail in 

MM 



We also note that the observed velocity dispersions of 
the most extreme outliers from the size-mass relation 
(Alogi? c < —0.4) are only marginally higher than for 
galaxies with 'normal' sizes. Only one of these candidates 
(logM* = 10.73) has A log a > 0.18 dex; the median 
value of A log a is 0.03 dex. It would seem that the effects 
of 'outlier noise' (ie. objects being pushed to the edge of 
the observed distribution by measurement errors, rather 
than their true, intrinsic properties) become dominant at 
these very extreme values of A log R . 

With these caveats, the observed velocity dispersions 
generally support the idea that the offsets from both the 
M* — R c and M* — a relations for our compact galaxy 
candidates can be explained by their having small sizes 
for their masses/velocity dispersions. 

4. THE IMPORTANCE OF SELECTION EFFECTS 
FOR COMPACT GALAXIES 

4.1. SDSS Spectroscopic Sample Selection 

In order to be targeted for SDSS spectroscopic follow- 
up (and so to appear in Figure [2]), galaxie s have to sat- 
isfy a complicated set of selection criteria ([Strauss et al.l 
2002). In brief, there is a magnitude cut: objects must 
be detected at > 5a significance, and have rp < 17.77. 
Any objects that have been marked as blended and then 
segmented into smaller objects are rejected, as are any 
objects that include saturated pixels, or have been de- 
blended from objects with saturated pixels. There are 
also a series of (low) surface-brightness-dependent crite- 
ria that are not relevant here. 

The first important selection criterion for our purposes 
is the star/galaxy separation criteria, since we are con- 
cerned about bright, compact galaxies being mistakenly 
identified as stars. Star/galaxy separation is done on 
the basis of the difference between the 'PSF' and 'model' 
magnitudes in the r band. (Here, the PSF magnitude 
is derived by fitting the PSF model to each object, in 
analogy to the exponential/de Vaucouleurs model fits 
described in §2.1[ and then aperture corrected to 7'.'4.) 
Specifically, objects are only selected for spectroscopic 
follow-up where: 

A SG = r PSF - r mode i > 0.3 . (2) 

Further to this star/galaxy discriminator, in order to 
avoid cross talk between spectroscopic fibers, galaxies 
with fiber magnitudes g < 15, r < 15, and i < 14.5 are 
also rejected. Lastly, all objects with rp < 15 and a Pet- 
rosian radius Op < 2" are rejected. This criterion was 
introduced to eliminate "a small number of bright stars 
that that managed to satisfy equation [2] during the com- 
missioning phase of the survey, when the star/galaxy sep- 
aration threshold was Age — 0.15 mag, and was retained 
for later runs to avo id saturating the spect roscopic CCDs 
([Strauss et alj|2002| ). IStrauss et all (|2002t ) also note that 
of the approximately 240000 r < 17.77 objects in runs 
752 and 756, none were rejected by the rp < 15, Op < 2" 
criterion alone. 

In order to model these selections, we need to relate 
the relevant observed quantities (viz., the apparent Pet- 
rosian magnitude, rp, gri fiber magnitudes, the apparent 
Petrosian size, Op, and the star/galaxy separation pa- 
rameter, A§g) to intrinsic size and stellar mass. 

For a given redshift/distance, the intrinsic size can be 
trivially related to the observed effective radius, O c . In 
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order to relate rp to M* , we have made a simple fit to the 
relation between stellar mass and absolute magnitude in 
the observers frame r band [ie. with no K-correction) 
for red sequence galaxies at 0.066 < z < 0.12 with 
M r > —21. Note that this method naturally accounts 
for mass-dependent trends in, e.g., metallicity along the 
red sequence. The scatter around this relation is ~ 0.06 
dex, with no obvious magnitude dependence. We have 
derived similar relations for both gp and ip. 

We have also derived empirical relations for 8p, Asg, 
and the difference between the Petrosian and fiber mag- 
nitudes, Afib = mp — mfibor, as functions of © e and rp, 
using the sample of massive, red sequence galaxies shown 
in Figure^ The scatter around these relations is 0.059 
dex (15 %) , 0.18 mag (18 %), and 0.11 mag (9 %) re- 
spectively, with no obvious systematic residuals. 

Note that there is a danger of circularity in this argu- 
ment: any objects that do not satisfy the selection cri- 
teria will not be present in the sample that we are using 
to model the selection criteria. The crucial assumption 
here, then, is that we can extrapolate the functions for 
Op (So, r P ), A SG (O c , r P ), and A fib (©c, r P ) down past 
the limits of the spectroscopic sample. In this regard, it 
is significant both that the derived functions are smooth 
all the way down to the selection limits, and that we do 
not see obvious cut-offs in the data associated with these 
limits. 

In Figure [21 we show how these selection criteria trans- 
late onto the (M*,R e ) plane for several example redshifts 
between 0.035 and 0.12. The thicker, roughly diagonal, 
long-dashed lines show the star/galaxy separation crite- 
rion; the dotted lines show the 'cross-talk' fiber magni- 
tude selection; the thinner, short-dashed boxes show the 
effect of the 'saturation' selection against bright, compact 
objects. Note that, for example, a galaxy with M* > 10 11 
M Q and R c < 2 kpc would not be selected as an SDSS 
spectroscopic target for z < 0.05. 

4.2. Compact Galaxies at High and Low Redshifts 

In Figure [6] we again show the size-mass relation for 
our sample of massive, red sequence galaxies at 0.066 < 
z < 0.12, with the exception that we have not plotted 
those galaxies rejected as per ij3.ll Furthermore, in con- 
trast to Figure [2j we have used the selection limits de- 
rived in £!4.1l to estimate the relative completeness of the 
SDSS spectroscopic sample across the 0.066 < z < 0.12 
volume; these are shown by the contours. These com- 
pleteness estimates also include the rp < 17.77 selec- 
tion limit, which can be seen to affect galaxies with 
M* < 10 10 6 Mq at the distant end of our redshift win- 
dow. 

For comparison, we h ave also overplotte d the high- 
redshift samples of ID09I (yellow points) and IvDOSl blue 
(blue points). Where we have used size measurements 
from the z-band for the SDSS galaxies, these high- 
redshift studies use the NICMOS F160W filter, which 
corresponds to rest- frame r at z = 1.6, moving close to 
g by z — 2.3. Locally, the difference between z- and r- 
band measured sizes leads to a slightly different slope to 
the size-mass relation for red sequence galaxies (a slope 
of 0.65, rather than 0.56). The r- and z- band size-mass 
relations intersect at around M* ~ 10 10 M ; the mean 
r-band size at M* ~ 10 11 Mq is 15 % larger than in the 
z— band. That is, by using z-band derived effective radii, 



we are, if anything, underestimating the sizes of the local 
galaxies in comparison to those at high redshift. Sim- 
ilarly, our decision to use the De Vaucouleurs effective 
radii given in the basic SDSS catalog, rather than more 
general Sersic ones appears to lead to an underestimate 
of galaxy sizes. In other words, adopting r- or g-band 
derived sizes, or using Sersic instead of De Vaucouleurs 
effective radii, would increase the discrepancy between 
the high- and low-redshift samples. 

There is one of our candidates (marked with a cross) 
that appears to ha ve similar properties to one of the 
larger of the lvD08l galaxies. This turns out to be the 
galaxy whose observed velocity dispersion is inconsistent 
with its being genuinely compact ( §3.3|) : where we would 
predict A log a — 0.24 dex, or a = 310 ± 70 km/s, what 
we observe is A logo- = —0.17 dex and a = 129 ± 14 
km/s. This is also the galaxy with the largest dif- 
ference between the Sersic- and De Vaucouleurs-sizes 
(log Rucv/Rscr = -0.34; see %T2|) . Adopting the NYU 
VAGC Sersic size measurement is not sufficient to rec- 
oncile the observed size and mass with the velocity dis- 
persion: the observed velocity dispersion would still be 
too small by ~ 0.2 dex. This galaxy also sits nearly 0.25 
dex above the median c olor-mass-to- lig h t relat ion shown 
in Figure [TJ using the iBell fc de .Tonal (f200lh prescrip- 
tion for M*/L as a function of (B — V) leads to a stellar 
mass estimate that is 0.17 dex lower. Adopting both this 
mass estimate and the NYU VAGC size estimate, we do 
find consistency between Alogi? e and Aloger. In this 
sense, this galaxy is the weakest of our compact galaxy 
candidates — it seems to have had its size underestimated 
and/or its mass overestimated. 

We also stress that the observed velocity dispersions of 
the candidates that lie furthest from the main size-mass 
relation suggest that these galaxies have had their sizes 
significantly underestimated (see $3.3$ . 

If the lvD08l galaxies were placed at 0.066 < z < 0.12, 
the SDSS spectroscopic completeness would typically be 
> 75 %. Note, ho wever, that there are two R e < 0.5 kpc 
galaxies from the lvD08l sample and one from the [D09 
sample for which the SDSS completeness is jus t 20-40 
%. The average SDSS completeness for the lvD08l galaxies 
placed a t 0.066 < z < 0.12 would b e 80 %. 

If the iKriek et alJ (j2008aj) / fvD08l galaxies were not to 
evolve in either size or number density from z ~ 2 to 
the present day, we would expect there to be ~ 6500 
M» > 10 11 M galaxies with Alogi? c < -0.4 dex at 
0.066 < z < 0.12, of which ~ 5250 should appear in the 
SDSS spectroscopic sample. Instead, we have only one 
weak candidate. 

As an interim conclusion, then, we have shown that 
there are no galaxies in the local universe (at least as 
probed by the SDSS spectroscopic sample) that are di- 
rectly analogous to the compact galaxies found at high 
redshift. This dearth of compact galaxies cannot be ex- 
plained by selection effects. In Appendix [Al we confirm 
this conclusion by searching for compact galaxy candi- 
dates from within the SDSS photometric sample, using 
photometric redshifts. 

Moreover, we stress that those galaxies which we have 
identified as 'compact' are not qualitatively similar to 
the compact galaxies found at higher redshifts, which are 
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Fig. 6. — The size-mass relation for massive, red sequence galaxies at low and high redshifts. The points and circles are for SDSS 
galaxies, as in Figure[2] those galaxies that we have rejected as per ij3]are not shown. The contours show the relati ve volume completeness 
of the SDSS spectroscopic sample for 0.066 < z < 0.12, as marked. The orange points with error bars are the D09 sample of 1.2 < z < 2.0 
galaxies from the GDDS and MUNICS. The blue points with error bars are the lvDOa sample of z ~ 2.3 galaxies from MUSYC. There are 
no red galaxies in the local universe with sizes and masses comparable to the compact galaxies found at high redshift. This lack cannot be 
explained by selection effects. 



offset from the local size-mass relation by at least twice 
as much again as our local compact galaxy candidates. 

4.3. The Number Density of Massive, Compact 
Galaxies 

In Figure [7J we provide a more quantitative statement 
of our conclusion with respect to the size evolution of 
massive galaxies from z ~ 2 to z ~ 0.1 by plotting 
the size distribution for massive, red galaxies in different 
mass bins. In this figure, the filled histograms represent 
the main SDSS spectroscopic sample described above. 
The horizontal-hatched histograms show, for compari- 
son, the situation at z ~ 1.6, based on the ten ID09I 
galaxies drawn from the GDDS; similarly, t he diagonal- 
hatched histograms show the n ine z ~ 2.3 iKriek et al.l 
(|2008af) galaxies wi th sizes from lvD0"8l 

The (Kriek et alJ (|2008af ),/ ivD0i sample is representa- 
tive, but not complete. In order to derive the de nsities 
plotted in Figure [7J we have scaled each of the lvD08l 



galaxies as follows: first, we have normalized the distri- 
bution to have a density of 1.5 x 10~ 4 Mpc~ 3 , which 
corresponds to the total nu mber density of all 2 < z < 3 
galaxies to the mass limit of lKriek et al.l (l2008al). derived 
using the mass function fit given by iMarchesini et afl 
(2008); then, we have scaled this distribution by a fac- 
tor of 16/36 to count only tho se galaxies with lit tle or 
no ongoing star formation from lKriek et al. (2008a[) that 
see m to form a red sequence (IKriek et al.ll2008bfl . For 
thelD09l sample, we are able to use 1/V ms , x scalings from 
iGlazebrook et alJ 1)20040 . 

The location of each individual high-redshift galaxy is 
marked in Figurc[7jwith an arrow: the slightly lower blue 
arrows show the[vp08 galaxies; the slightly higher yellow 
arrows are for the[D09 galaxies. Clearly, given the small 
numbers, the uncertainties on these high redshift values 
are quite large, but they do provide a useful order of 
magnitude estimate for comparison to the local values. 
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Fig. 7. — The observed size distribution of massive, red galaxies 
at z ~ 0.1, Z ~ 1.6, and z ~ 2.3; each panel is for a different mass 
range as marked. In each panel, the solid histogram represents 
the SDSS spectroscop ic sam ple. The blue, diagonally-hatched his- 
togram represents the vD08 sample of nine massive, passive galax- 
ies at z ~ 2.3; the yellow, horizontally-h atche d histogram repre- 
sents the 10 z ~ 1.6 GDDS galaxies from lD0"9l . The arrows at the 
bottom of each panel indicate the positions of the individual high 
redshift galaxies. The z ~ 2.3 galaxies have to undergo signifi- 
cant structural evolution over z < 2.3 to match the properties of 
local universe galaxies; at least part of this evolution has already 
occurred by z ~ 1.6. 



The clear implication from the comparison between the 
z ~ 0.1 and z ~ 2.3 d ata in Figure [3 is tha t, con sistent 
with the conclusions of lvD08l , not one of the|vl)08 galax- 
ies is consistent with the properties of the z ~ galaxy 
population. With the results we have now presented, we 
can extend this conclusion by confirming that this dis- 
crepancy cannot be explained by selection effects in the 
low redshift sample. 

There are local analogs for less than half of the z ~ 1.6 
galaxies, albeit with considerably higher number densi- 
ties. This would imply that at least some (< 50 %) of 
the z < 2.3 evolution has already occurred by z ~ 1.6. 

5. DISCUSSION 
5.1. Compact Galaxy Properties 

In Figure we plot the properties of our compact 
galaxies in comparison to the general massive, red galaxy 
population. In each panel, the circles highlight our com- 
pact galaxies, while the points and greyscale show all 
galaxies with M* > IO 10 7 M , 01 (u - r) > 2.5, and 
0.066 < z < 0.12. The large grey boxes with error bars 
show the mean and standard deviation of each plotted 
property in quintiles of the velocity dispersion distribu- 
tion. Similarly, the red boxes with error bars show the 
mean and standard deviations for our compact galaxy 
candidates in two bins, separated at a = 200 km/s; the 
median for this sample. 

In each of the panels of Figure[8](from left to right), we 
show the equivalent width of the US line (where negative 
values imply emission), the luminosity weighted mean 
stellar metallicity, and the luminosity weighted mean 
stellar age, as deriv ed from the DR4 SDSS spectra by 
iGallazzi et al.l (|2005h . Because these estimates are avail- 
able only for DR4, only around half of our compact can- 
didates can be plotted in these panels, and only 3/10 of 
those with M* > 10 11 M . 

We have also matched our compact galaxy s a mple t o 
the AGN sample described bvlKauffman n~et al.l (|2003bl) . 
for SDSS DR4. These AGN hosts have been selected by 
their [OIII1/H/3 and [Nil] /Ha e missio n line ratios; ie. the 
iBaldwin. Phillips fc Terlevichl (|1981l . BPT) diagram. 34 
of our 63 galaxies appear in the DR4 catalog; of these, 
11 are classified as AGN on the basis of their emission 
line ratios. This is slightly higher than the AGN fraction 
of the parent sample, which is in the range 20 — 26 % for 
the mass range we are considering. Of the 11 galaxies 
identified as AGN hosts, five sit on or slightly above the 
main M* — £[oin] relation, with Lroini ~ 10 6 L Q , four 
have i[oni] ~ 10 7 ~ 8 L Q , and one is quite high luminosity, 
with L[oin] = 10 87 Lq. These 11 galaxies are marked in 
e ach panel of Figur e jH with a small blue cross. 

iKauffmann et alj (|2003bf) also provide revised stellar 
mass and velocity dispersion measurements for these 
galaxies. Accounting for the presence of an AGN does 
not have a major impact on these measurements: the 
masses and velocity dispersions change at the level of 
0.05 dex and 16 km/s, respectively. That is, while it is 
possible that an optically bright point source may bias 
the measured sizes of these galaxies downwards, within 
the stated errors, the AGN does not significantly affect 
the derived values of M* or a. (It is relevant here that 
only one of our compact galaxy candidates shows a signif- 
icant residual point-source after subtracting off the best- 




Fig. 8. — The properties of our compact galaxy candidates in comparison to the general population of massive, red sequence galaxies. 
As in other Figures, the black points and greyscale show the data density of all galaxies with M„ > 10 10 7 Mq, 01 (m — r) > 2.5, and 
0.066 < z < 0.12; the red circles show our compact galaxy candidates. The grey boxes with error bars show the mean and rms scatter in 
each quantity for quintiles of the velocity dispersion distribution; the red boxes show the same for those of our compact galaxy candidates 
with a < 200 km/s and a > 200 km/s separately. At fixed velocity dispersion, our compact galaxies have slightly lower than average mean 
ages and slightly higher metallicities — however, this result is only significant at the 2a level. 



fit Sersic profile, as produced by GALFIT; see j)3.2l ) 

Looking now at Figure [SJ it is clear that the majority 
of our compact galaxy candidates have quite old stellar 
populations. For the a < 200 km/s bin, the median 
age is 6 Gyr, although the ages do range from 2 to 10 
Gyr, while all but one of the a > 200 km/s candidates 
have T > 6 Gyr. Among the lower velocity dispersion 
candidates, there is a clear tendency towards relatively 
high equivalent widths for H<5 absorption, suggestive of 
a relatively recent (< 2 Gyr) star formation event. 

At fixed velocity dispersion, our compact galaxy 
candidates may have slightly higher metallicities, and 
be slightly younger than average. Using bootstrap- 
resampling on similar sized samples drawn randomly 
from the mass-limited sample, and controlling for ve- 
locity dispersion, these results are weakly significant at 
best: 1.9a for the age, and < la for metallicity. Con- 
sidering only the higher velocity dispersion candidates 
(a > 200 km/s), the significance of these differences be- 
come 1.9a and 2.2u for the age and metallicity offsets, re- 
spectively. This weakly sign ificant result should b e con - 
tras ted with the resul t s of Shanka r fc Bernardil (2009) 
and Ivan der Wei et alj (|2009fl . who find that, on aver- 
age and at fixed dynamical mass, early type galaxies 
with higher velocity dispersions (or, equivalently, smaller 
sizes) have older mean stellar ages. 

While the younger mean stellar ages and lower metal- 
licities of our compact galaxy sample are only weakly 
significant, both would imply a relatively late start to 
star formation for these galaxies and/or their progeni- 
tors. But if these galaxies grow in size through mergers 
(for example) then it is possible that these galaxies are 
small not because their formation is delayed relative to 
other galaxies of the same mass or velocity dispersion, 
but rather because they have had fewer mergers overall, 
or perhaps just fewer recent mergers. That is, it may be 
that, at fixed mass, these compact galaxies are in fact 



older, in the sense that they have been assembled ear- 
lier, and existed in (more or less) their present form for 
longer than other galaxies of the same mass or velocity 
dispersion. 

5.2. Comparison to Other Recent Works 

In a similar study to this, using sizes and photome- 
try fro m the NYU VAGC for SDSS DR6, ITruiillo et al.1 
(2009) have recently reported the detection of 29 z < 0.2 
galaxies with M, > 8 x 10 10 M© and R c < 1.5 kpc. 
In contrast, we find just one galaxy from our 0.066 < 
z < 0.12 red sequence galaxy sample that satisfy these 
mass and size criteria; this implies a difference in volume 
densities of a factor of 5.5. Most of this difference is ex- 
plained by the fact that we have preselec t ed ou r compact 
galaxies to be red. Of the ITruiillo et al.l (|2009h galaxies, 
only around 30 % (9/29) satisfy our 01 (u - r) > 2.5 
criterion, bringing our number densities into agreement. 
On the other hand, if we also look at 01 (u — r) < 2.5 
galaxies, inspected as per ^3.11 we find only 7 additional 
candidates, only 3 of which have M» > 8 x 10 10 M Q ; the 
most massive of these blue compact galaxy candidates is 
7.2 x 10 10 M Q . 

We also note that the ITruiillo et al.l (|2009f ) galaxies 
have considerably smaller observed sizes than the galax- 
ies we consider here. (Recall that we adopt a minimum 
observed size of 0'.'75; for galaxies with inferred sizes 
smaller than this, we adopt a size of 0'.'75 as an upper 
limi t on the true s ize.) T he largest observed size among 
the ITruiillo et all (|2009f ) sample is 0'.'70; the median is 
just 0'.'48. Enf orcing our minimum allowed size of C/75, 
only one of the lTruiillo et ahl (|2009f ) galaxies, irrespective 
of color, would have R c < 1.5 kpc; the median size for 
the sample would bec ome 2.1 kpc. I t is als o relevant here 
that only one of the ITruiillo et alj (|2009f ) galaxies is at 
z < 0.12; the other 28 are all found at 0.1 2 < z < 0.20. 
This suggests that the lTruiillo et all (|2009h size measure- 
ments may be biased by inadequate resolution. 
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The key diffe rence between our co mpact galaxy sam- 
ple and that of ITruiillo et al.l ([2009) is that they find a 
median velocity dispersion which is only 0.04 dex higher 
than their control sample, eve n though the mea n size 
and mass are offset from the IShen et al.1 (|2003f ) rela- 
tion by —0.5 dex. This discrepancy can only be ex- 
plai ned by either very l arge structural differences, or if 
the ITruiillo et alj (|2Q09h sample is disproportionally af- 
fected by large measurement errors in size and/or mass. 
In this context it is significant that, the observed velocity 
dispersions for our candidates with Alogi? e < —0.4 dex 
imply that their very small inferred sizes are produced by 
large errors in the measured sizes (see §3.3[) . Only with 
followup observati ons will we be able to determine the 
true nature of the I Truiillo et ail (|2009f ) galaxies. In any 
case, ITruiillo et al.l (|2009f ) also do not find any galaxies 
directly comparable to those found at z > 2. 

Even more recently, IValcntinuz zi et ail (|2009T ) have de- 
scribed a sample of 147 compact galaxies selected from 
the Wide-field Nearby Galaxy-cluster Survey (WINGS) 
of X-ray sele cted clusters at 0.04 < z < 0.07. Unlike 
in this work, IValcntinu zzi et aT] (|2009f) do find g alaxies 
with properties comparable to the 3 largest |vD08 galax- 
ies; similarly, there are local WINGS analogs for 8 of 
the 10 GDDS galaxies from ID09I . However, this relies 
on their scaling the high redshift galaxies' masses down 
by 0.15 dex to account for the poor treatment of NIR- 
luminous thermally p ulsating asymptotic giant branch 



(TPAGB) stars in the iBruzual & Charlotl (1200 3) stellar 
population models. While both iKriek et al.l (|2009f ) and 
M uzzin et all (|2009fl show that the stellar masses for the 
z ~ 2.3 implied by different models vary by ~ 0.1 dex, 
we have not applied such a correction here. 

We note, however, that the high- and low-redshift sam- 
ples have been treated consistently here, including the 
fact that all masses were deriv ed from the rest-frame op - 
tical. Further, we note that Ivan der Wei et all (12009 ) 
have shown that stella r masses derived fr o m th e rest- 
frame optical and using IBruzual fc Charlotl (|2003l ) mod- 
els are consistent with the dynamical masses of z < 1 
galaxies, and are unaffected by the TPAGB uncertain- 
ties. 

The IValentinuzzi et all (|2009h compact galaxies sam- 
ple is selected by effective surface mass density, S = 
M*/2-kRI > 4x 10 9 Mq kpc~ 2 , in the range 3—50 xlO 10 
M Q . Our compact galaxy selection is roughly equivalent 
to S c > 3.6 x 10 9 M Q kpc" 2 ; nearly half (28/63) of our 



compa ct galaxy candidates satisfy the Valentinuzzi et al 
2009) S c criterion. For their sample. [Valentinuzzi et al 



20091 ) derive a number density of 1.6 x 10~ 5 Mpc -3 ; to 
our mass limit of M* > 10 107 M Q , this value becomes 
1.2 x 10 -5 Mpc~ 3 . These values are solid lower limits, 
as they assume that no such galaxies exist outside of the 
clusters observed by WINGS. For our sample, however, 
the number density of M* > 10 107 M Q galaxies with 
Eo > 4 x 10 9 M kpc" 3 is just 2.85 x 10" 7 Mpc' 3 . 

That is, after correcting as best we can for the differ- 
ent stellar mass limits, our number densities are incon- 
sistent by a factor of mo re than 40 with those found by 
IValentinuzzi et all (|2009h . Again, our use of z-band effec- 
tive radii leads to smaller measured sizes than for bluer 
bands; this discrepancy would only increase using r- or 
g-band measured sizes. Either our results are badly af- 



fected by unexplained selection effects, or there are large 
discrepa ncies between our size an d mass estimates and 
those of lValentinuzzi et"aH (|2009h . 

We have considered possible spectroscopic selection ef- 
fects that could bias against bright, compact objects 
in §4.11 and shown these to be rele yant for z 5, 0.05. 
These effects may well explain why IValentinuzzi et ail 
(2009) were able to match only a small fraction of their 
compact galaxies (which have 0.04 < z < 0.07) to ob- 
jects in the (DR4) SDSS spectroscopic catalog. We have 
shown, however, that our 0.065 < z < 0.12 results 
are not strongly affected by these kinds of selection ef- 
fects (Figure [6j see also Appendix A). Our estimated 
completeness is more th an 60 % for all galaxies in the 
IValentinuzzi et ail (|2009h sample and greater than 90 % 
for 90 % of the sample. The selection effects considered 
in £14.11 thus cannot explain the difference in our inferred 
number densities. 

An alternativ e ex planation is that the 
IValentinuzzi et ail (|2009l ) galaxies only exist in rich 
clusters, and that SDSS suffers much higher spectro- 
scopic incompleteness in such dense fields because of 
fiber collisions. A comple tely indepdendent estimate 
can be obtained from the iFaber et al.l (|1989l ) sample: 
we find that 5/319 of these galaxies have sizes smaller 
than the mass-size relation by a factor of 2 or more. 
This fraction for clusters that is approximately 15 times 
higher than what we find for all galaxies in SDSS. 
While not conclusive, this does suggest that SDSS 
may suffer from additional incompleteness beyond the 
effects we con sider here. That said, w e note that several 
studies (e.g. iKauffmann et al.l l2004t iPark et alj 120071 : 
iWeinmann et al.l I2009T ) have found little or no evidence 
for an environmental dependence of the size-mass 
relation within SDSS. 

Thus we can find no easy explanati o n for the dif- 
ference between the IValentinuzzi et all (|2009l ) results 
and our own. Here again, velocity dispersion measure- 
ments would provide an useful consistency check on the 
IValentinuzzi et al.l (2009) size and mass measurements. 

Even despite these diff erences, however, we note that 
IValentinuzzi et ail (|2009f ) conclude that — barring large 
systematic errors in the high-redshift measurements — at 
least 65 % (cf. our value of 100%) of the z ~ 2.3 galaxies 
from lvD08l and at least 20 % (c f. our value of 60 %) of 
the z ~ 1.6 galaxies from [P09l have disappeared from 
the local universe. Accepting the high-redshift results, 
these galaxies simply cannot evolve passively and stat- 
ically into the red sequence and/or early type galaxies 
found in the local universe. 

6. SUMMARY AND CONCLUSIONS 

The central question of this work has been the exis- 
tence or otherwise of massive, compact, quiescent and/or 
early type galaxies in the local universe, and particularly 
the importance of selection effects in the SDSS spectro- 
scopic sample for such galaxies. We have shown that, 
especially for lower redshifts (z < 0.05), galaxies with 
the masses and sizes of those found at z > 2 would not 
be targeted for spectroscopic followup (Figure H]) . The 
main reason for this is not the star/galaxy separation 
criterion, but rather the exclusion of bright and compact 
targets in order to avoid saturation and cross-talk in the 
spectrograph (see ^4. II) . 



14 



Taylor, Franx, Glazebrook, Brinchmann, van der Wei & van Dokkum 



We have therefore conducted a search for massive, com- 
pact galaxies at 0.066 < z < 0.12, where these selection 
effects should be less important. We estimate that for 
0.066 < z < 0. 12, th e av erage completeness for galaxies 
like those from lvD08l and lD09l would be > 25 % at worst, 
and ~ 80 % on average (Figure [6|) . 

Starting from a sample of massive (M* > 10 10 7 M©) 
red sequence ( 01 (w — r) > 2.5) galaxies, we have se- 
lected the 280 galaxies with inferred sizes that are a fac- 
tor of 2 or more s maller than would be expected from the 
IShen et all (|2003[ ) M*-R e relation for early type galaxies. 
In order to confirm their photometry and size measure- 
ments, we have visually inspected all of these objects. 
Unsurprisingly, by selecting the most extreme outliers, a 
large fraction of these objects (~ 70%) appear to be in- 
stances where the size and/or stellar mass estimates are 
unreliable f ij3.ip . 

For the 63 galaxies with no obvious reason to suspect 
their size or stellar mass estimates, there is good agree- 
ment between the default SDSS size measurement (based 
on the 2D light distribution, using a sector-fitting algo- 
rithm, and assuming a de Vaucouleurs profile), and those 
given in the NYU VAGC (based on the azimuthally av- 
erage growth curve, assuming a more general Sersic pro- 
file). However, particularly for galaxies with high n, 
the de Vaucouleurs size measurement is systematically 
smaller than the Sersic one, at the level of < 25 % (j j3.2j) . 

In general, and as expected, our 63 compact galaxy 
candidates have significantly higher than average veloc- 
ity dispersions (Figured]). While it remains possible that 
the sizes of at least some of our compact galaxy candi- 
dates may have had their sizes underestimated by ~ 30 
%, in general, the relatively high observed velocity dis- 
persions support the notion that they are indeed unusu- 
ally compact given their stellar masses ( 33, 3[) . 

Among our compact galaxy candidates, there are no 
gal axies with sizes comparable to those fou nd z ~ 2.3 
by lvD08l ; we find analogs for < 50 % of the ID09I galax- 
ies at z ~ 1.6 (Figures [6] and EJ. This lack cannot be 
explained by selection effects. To confirm this, we have 
also compared the size-mass diagram, constructed using 
photometric redshifts, based on both the full photomet- 
ric sample and the spectroscopic sub-sample (Appendix 
[A"l) . While it is conceivable that SDSS is missing a few 
massive, compact galaxies, the re are a gain n o signs of 
galaxies comparable to those of vD08l or lD09L 

It is not impossible that some systematic errors in the 
estimation of M*/Ls for the high redshift galaxies (e.g., 
an evolving stellar IMF) mean that their stellar masses 
are vastly overestimated, however it wo uld req uire an 
overestimate of > 0.7 dex to reconcile the lvD08l galaxies 
with the sizes of the smallest galaxies we have identified 



in the SDSS catalog. 

Accepting the high redshift observations at face value, 
then, our results confirm that massive galaxies, both indi- 
vidually and as a population, must undergo considerable 
structural evolution over the interval z < 2.3 in order to 
develop into the kinds of galaxies seen locally — even after 
star formation in these galaxies has effectively ended. We 
see some hints that a significant amount of this evolution 
(< 50 %) may have already occurred by z ~ 1.6. 

The fact that each and every one of the lvD08l galaxies 
must undergo significant structural evolution to match 
the properties of present-day galaxies implies that the 
mechanism that drives this growth must apply more or 
less evenly to all galaxies. To see this, let us assume 
that some external process drives the size evolution of 
these galaxies, and that even a single event is sufficient 
to move an individual galaxy onto the main size-mass 
relation. Then, we can assume some simple probabil- 
ity distribution for the number of events, N, among in- 
dividual galaxies. (For example, we could assume that 
events occur randomly across the time interval z < 2.3, 
or that each galaxy experiences N ± VA eve nts.) N ow, 
our results suggest that the number density of lvD08l -like 
galaxies drops by at least a factor of 5000 since z ~ 2.3. 
In order to ensure that at most 1/5000 galaxies have 
N = after z ~ 2.3, simple probabilistic arguments im- 
ply that the average galaxy must undergo > 20 events. 
This would imply that a strongly stochastic process like 
major mergers cannot be the primary mechanism for the 
strong size evolution of massive galaxies. 

Apart from their small sizes and high velocity disper- 
sions, our compact galaxy candidates are not obviously 
distinct from the general population (Figure [H]). If any- 
thing, at fixed velocity dispersion, our compact galaxies 
have stellar populations that are slightly younger than 
average (at ~ 1.9er significance). Even so, the majority 
of these galaxies' stellar populations are definitely 'old', 
with luminosity-weighted mean stellar ages typically in 
the range 6-10 Gyr. But if some external mechanism 
drives the size evolution of these galaxies, we speculate 
that their small sizes may indicate that they have as- 
sumed their present form comparatively early, and in 
this sense they may actually be relatively old (see also, 
e.g.. Ivan der Wei et alj '2009). If so, with better under- 
standing of the processes that determine the sizes of early 
type galaxies, and in particular the role of merging, the 
properties of these galaxies could provide a means of con- 
straining the evolution of massive galaxies after they have 
completed their star formation, including their late-time 
merger histories. 



APPENDIX 

A. LOOKING FOR MASSIVE COMPACT GALAXIES IN THE SDSS PHOTOMETRIC SAMPLE 

In this Appendix, we present a complementary analysis in which we directly compare the spectroscopic and photo- 
metric samples, in order to test the conclusion that the lack of massive, compact galaxies in the spectroscopic sample 
cannot be explained by the selection effects. 

A.l. Selecting Galaxies by Color Alone 

Before we can address the question of massive compact galaxies in the SDSS photometric sample, we must first 
devise a means of separating stars and galaxies without selecting on the basis of observed size or light profile. Our 
method for doing so is shown in Figure [H which plots the observed (extinction-corrected) ugrz colors of different 
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Fig. 9. — Selecting z < 0.1 galaxies based on color alone. Each panel shows the observed (u — g)—(r — z) color— color plot for objects in 
the SDSS spectroscopic sample; the right panel simply shows the central region in greater detail. Points are color-coded according to their 
spectral classification, viz.: galaxies [grey), galaxies with 0.066 < z < 0.10 (light red), galaxies with 0.066 < z < 0.10 and 0,1 (u — r) > 2.5 
(dark red), quasars (yellow), late-type stars (light blue), and ordinary stars (dark blue), . The box shows the color selection that we use 
to select z < 0.1 galaxies. This selection should produce a reasonably complete sample of z < 0.1 galaxies, with some contamination from 
both stars and quasars. In particular, later-type stars and some quasars have observed SEDs that are very similar to red sequence galaxies 
at z ~ 0.1. 

classes of objects from the spectroscopic sample; we show: all galaxies (grey), 0.066 < z < 0.12 galaxies (bright red), 
and those with (u — r) > 2.5 (dark red), O-K stars (dark blue), M-type or later stars (light blue), and quasars 
(yellow) . 

The black box shown in Figure [9] shows our criteria for selecting 0.066 < z < 0.12 galaxies based on their ugrz 
colors: 

0.6 <{u- g) < 2.4 and 
0.3 x (u-g)<(r-z) < 1.2 . (Al) 

Again, we apply this selection in terms of model colors. Note how, whereas the stellar sequence is reasonably well 
separated from the region of color space occupied by galaxies for (u — g) < 2.5, beyond this point, the late- type stellar 
sequence turns up, such that late-type stars and galaxies are blended. In the most general terms possible, the mean 
galaxy redshift increases towards redder (u — g) colors. This means that our ability to distinguish red galaxies from 
late- type stars on the basis of their optical SEDs is limited to z < 0.12. 

In the right-hand panel of Figure [51 we zoom in on this selection region. From this panel, it is clear that a large 
proportion of quasars will also be included in our color-selected 'galaxy' sample. Similarly, it is clear that this color 
selection is not 100 % efficient in excluding stars from our sample: more quantitatively, with this selection we are able 
to exclude more than 80 % of spectroscopically identified stars that are given 0.066 < z p hot < 0.12, while retaining 
more than 97 % of all 0.066 < z p hot < 0.12 galaxies. Furthermore, it should be remembered that stars are already 
heavily selected against for the spectroscopic sample plotted in Figure [9l the relative number of stellar 'contaminants' 
may well be considerably higher for the photometric sample. 

A. 2. Photometric Redshifts and Stellar Mass Estimates 

A major improveme nt in DRY is a complete revision in how the basic (photoz) photometric redshifts are derived 
(|Abazaiian et aLll2009f) . Rather than using some combination of synthetic template spectra to reproduce the observed 
colors of individual galaxies, the new photoz algorithm directly compares the observed photometry of individual 
galaxies to that of galaxies that have spectroscopic redshifts. Specifically, for each individual object, the algorithm 
finds the 100 closest neighbours in ugriz color space, and fits a hyper-plane to these points, rejecting outliers; the 
redshift is then determined by interpolating along this 4D surface. In comparison to the DR6 algorit hm, this reduces 
the R MS redshift error by more than 75 % ((Az) = 0.025), and significantly reduces systematic errors (|Abazaiian et al.l 
f2009h . 

For this analysis, rather than full SED-fit stellar mass estimates assuming the photometric redshifts, we will simply 
make use of the empirical relation between 01 (g — i) color and M*/L (Equation [1]). In this way, we are able to recover 
the z spcc -derived, SED-fit M*/Xs of the sample of galaxies shown in Figure Q] to 0.045 dex (la); including the effects 
of photometric redshift errors, k-corrections, and M*/L errors, the total (la) error in M* is 0.13 dex. This should be 
compared to the median formal error on the original SED-fit stellar mass estimates, 0.10 dex; that is, the errors on M* 
based on photometric masses (adding these two errors in quadrature) are only about 60 % greater than those based 
on spectroscopic redshifts. 
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Fig. 10. — The size-mass plot for massive, red sequence galaxies at 0.066 < z < 0.12, based on the spectroscopic and the full photometric 
SDSS catalogs. Each panel shows the sizes and masses of galaxies based on, from top to bottom, the spectroscopic sample using spectroscopic 
redshifts, the spectroscopic sample using photometric redshifts, and the photometric sample using photometric redshifts; in each case, only 
those objects inferred to have (u— r) > 2.5, and 0.066 < z < 0.10 are shown. In panel 3, many more objects with inferred sizes < 0.3 kpc 
can be seen; these are largely stars misclassified (in terms of their photometric redshifts) as galaxies. For Af» < 10 10,8 Mq and R e < 10 — 2 
kpc, comparison between panels 2 and 3 suggest that there may be a few additional galaxies in the photometric sample that do not appear 
in the spectroscopic sample. 



A. 3. The Size Distribution of Massive, Red Sequence Galaxies 

In Figure [TU1 we show three size-mass diagrams corresponding to, from top to bottom: (a.) the spectroscopic sample, 
analyzed using spectroscopic redshifts; (b.) the spectroscopic sample, analyzed using photometric redshifts; and (c.) 
the photometric sample, analyzed using photometric redshifts. In all three cases, the only selections applied to each 
sample are on photometric type (to exclude optical artifacts, etc., we require either a star or galaxy type classification) 
and ugrz color (to exclude stars); then, as in Figures and [H we are only showing those galaxies inferred to have 
0.066 < z < 0.12 and 0A (u — r) > 2.5. Again, objects with measured sizes smaller than 0'.'75 are shown as upper 
limits, assuming a size of 0'.'75. When comparing these three different analyses, the difference between (a.) and (6.) 
shows the effect of using spectroscopic versus photometric redshifts, and the difference between (6.) and (c.) shows the 
difference between the SDSS spectroscopic and photometric selection. That is, the comparison between (&.) and (c.) 
gives a direct indication of the level of incompleteness in the spectroscopic sample. 

Looking at panels (a.) and (6.), it is clear that the use of photometric redshifts produces a considerably greater 
scatter in the size-mass diagram, including a rather large number of galaxies with inferred stellar masses of 10 12 
Mq or greater. There is a clear excess of unresolved objects with inferred stellar masses greater than ~ 10 11 M Q 
in panel (b.) in comparison to panel (a.) However, we already know from section [3] that there are no objects in the 
spectroscopic sample with these sizes and masses — these objects cannot be genuine compact galaxies. Of the 34 with 
inferred M* > 10 11 Mq, 16 of these objects are spectrally identified as being stars, and one as a quasar at z = 0.102. 
Of the 17 spectrally confirmed galaxies, all have |z p hot — z spcc \ ^ 0.02. Of these, 15 have had their redshifts, and 
thus stellar masses, seriously overestimated; the other two are at z > 0.12, and so have had their intrinsic sizes 
underestimated. 

Turning now to the comparison between panels (b.) and (a), the first point to make is that the excess of unresolved 
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Fig. 11. — The observed size distribution of massive, red galaxies at 0.066 < z < 0.12; each panel is for a different mass range as marked. 
In each panel, the solid histogram represents the SDSS spectroscopic sample, analyzed using spectroscopic redshifts. The black and red 
histograms are the SDSS spectroscopic and photometric samples, respectively, analyzed using photometric redshifts. We have visually 
inspected all objects with inferred M* > 10 Mg and AR C < —0.5 dex; not one of these objects is a plausible massive, comapct galaxy 
candidate. The fact that the shape of the red histogram does not differ significantly from that of the black histogram for AR C > —0.5 dex 
indicates that the spectroscopic sample is not significantly biased against compact galaxies. 




sources is even more pronounced. We have matched all of these objects to the 2MASS point source catalog in order to 
investigate their NIR colors. 90 % of these objects fall in the stellar region of the (J — K)-K color-magnitude plot; 
similarly, 80 % fall in the stellar region of a (g — z)—(J — K) color-color plot. 

Further, we have visually inspected the 434 objects with inferred A/* > 10 11 M and with sizes smaller than the main 
M*-R e relation by 0.4 dex or more. Roughly 70 % of these objects are obviously stars: 133 come from crowded Galactic 
fields covered as part of SEGUE; 126 are double stars; 49 have clear diffraction spikes and/or are clearly saturated. 
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Another 12 objects have been cross-matched with the USNO-B star catalog (within 1"), and have measured proper 
motions of 1 — 4"/yr. 19 objects are the central point sources of very large spiral galaxies; most of these are also 
found in the ROSAT and/or FIRST catalogs. We also note that there are 17 very small disk or irregular galaxies 
with red point sources at or very near their centers. Most of these also have proper motion measurements from the 
USNO-B catalog, and several are spectrally identified as late type stars; it seems plausible that these galaxies simply 
have foreground stars coincidentally superposed very near their centers. 

In short, of the 434 objects from the full photometric sample that, on the basis of photometric redshifts, are inferred 
to have M* > 10 11 M Q and Ai? < —0.4 dex, not one remains as a viable compact galaxy candidate. 

A. 4. Estimating the Importance of Spectroscopic Selection Effects 

The conclusion from both the analyses that we have now presented is that there are no galaxies in the local universe 
with sizes and masses comparable to the compact galaxies found at higher redshifts. In Figure [TT1 we provide a more 
quantitative statement of this conclusion, by plotting the size distribution for massive, red galaxies in different mass 
bins. 

In this figure, the filled histograms represent the main SDSS spectroscopic sample, analyzed using spectroscopic 
redshifts, as in $3] The heavy black and red histograms represent the spectroscopic and photometric samples, respec- 
tively, analyzed using photometric redshifts, as in E1A.3I In all cases, objects excluded on the basis of visual inspection 
are not included; this accounts for the sharp cutoffs at A log R c = —0.3 and at A log R c — —0.4 for the filled and open 
histograms, respectively. Immediately above these cutoffs, where we have not visually inspected individual objects, 
but where there is likely to still be significant contamination, these distributions should be regarded as upper limits 
on the true distribution. In the upper panel, we plot those objects with observed sizes smaller than 0'.'75 separately 
as the light grey filled histogram, and the thin black and red histograms. 

As in the difference between the filled and solid black histogram, both of which are derived from the spectro- 

scopic sample, shows the increased scatter due to the use of photometric redshifts. 

Similarly, the difference between the black and red histograms show the difference between the spectroscopic and 
photometric samples, and so allow a quantification of the bias in the spectroscopic sample. By simply tallying the 
numbers of galaxies with —0.4 < Alogi? c < —0.3, we find that the 'completeness' (the ratio between the number of 
galaxies in the spectroscopic sample compared to the full photometric sample) is 75 %, 68 %, 67 %, and 43 % for each 
of these mass bins, from lowest to highest. 

In order to improve on these estimates, we have done the following. Using the approach described above, we have 
assigned each object a weight according to its z p hot - derived mass and size. Then, going back to the spectroscopic 
sample, we use these to compute the mean weight in cells of z spec -derived mass and size. The completeness contours we 
derive in this way are in good qualitative agreement with those shown in Figure[2j although they suggest incompleteness 
at the 2-5 % level for mean-sized galaxies with M* > 10 11 M Q . Using these values to estimate the number of M* > 10 11 
Mq galaxies with AR C < —0.3 dex, this suggests that the spectroscopic sample is missing on the order of 4 such galaxies. 
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